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EFFECTS OF AXIAL GRADIENTS OF VELOCITY AND MAGNETIC FIELD, 
ELECTRODE STAGGER, AND ION SLIP ON PARAMETERS 
IN MHD ACCELERATORS 

By Frank Hohl 
Langley Research Center 

SUMMARY 

The distributions of current and potential have been determined for two linear 
Faraday plasma accelerators that are under investigation at the Langley Research Center. 
The equations describing the electrical parameters of accelerators were reduced to two 
dimensions and numerical analysis was used to solve the equations. The internal resis- 
tance of the accelerator, the Hall potential gradient, and various other parameters have 
been calculated as functions of: velocity and magnetic field gradients along the channel 
axis, ion slip, various amounts of stagger of electrodes, the ratio of electrode width 
to insulator width, and the Hall parameter o> e r e . The results are compared with previ- 
ous theoretical work and with experimentally obtained data. 

INTRODUCTION 

In the Faraday accelerator, the Hall or electron drift current is reduced by using 
segmented electrodes which allow an electric Hall field along the channel axis. Since 
the electric field will be short circuited locally by the finite-length electrodes, strong 
nonuniformities in the current distribution can be expected near the electrodes. There- 
fore, a theory describing the electrical properties of an accelerator must use an analysis 
which is at least two dimensional. 

The mathematical formulation of the interaction between the electromagnetic fields 
and the plasma is given by the generalized Ohm's law which has been derived by Cowling 
(ref. 1). By using such an Ohm's law, the distribution of current and potential in plasma 
accelerators and generators has been derived by Hurwitz et al. (ref. 2), Schultz -Grunow 
and Denzel (ref. 3), Witalis (refs. 4 and 5), Dzung (refs. 6 and 7), and others by means of 
conformal mappings. Crown (ref. 8) and Celinsky and Fischer (ref. 9) have used numeri- 
cal methods to solve for the current flux function. Salvat (ref. 10) has used simple 
physical arguments to investigate the internal potential fall and resistance of magneto- 
hydrodynamic (MHD) generators. In most of these analyses, the effects of ion slip are 
neglected and in all the analyses, the fluid flow velocity and the applied magnetic field 



were assumed not to vary along the channel. Also the effect of shifting anodes with 
respect to the corresponding cathodes was not sufficiently investigated. 

In all accelerators there is a velocity gradient and in one of the accelerators that 
is being developed at the Langley Research Center (ref. 11) there is a relatively large 
gradient of the applied magnetic field in the flow direction. With a gradient of either 
velocity or magnetic field, the current and the electric field are no longer uniform in the 
channel center. The condition of homogeneous current distribution in the channel center 
which has been used by various authors (refs. 2, 4, and 5) who investigated the problem 
by means of conformal mappings no longer applies. Instead of uniformity in the channel 
center, some of the authors (refs. 3 and 12) assumed rectilinear current dividing lines 
oblique to the channel axis. This condition also is violated when the aforementioned 
magnetic field and velocity gradients are present. Of course, in addition to the non- 
homogeneity in the channel center, there are more serious problems if a solution is 
attempted by means of a conformal mapping in the presence of axial velocity and magnetic 
field gradients, since a Poisson equation must be solved instead of a Laplace equation. 

The purpose of the present analysis is to investigate the effect on the current and 
potential distribution of velocity and magnetic field gradients, of ion slip, of various 
amounts of stagger of the electrodes, of the ratio of electrode width to insulator width, 
and of the Hall parameter <u e r e . 

In the present analysis the boundary layers and the sheaths near the electrodes and 
insulators were neglected. Also the fluid properties are assumed not to vary across the 
channel. Kerrebrock (ref. 13) has investigated the effects of variations of the fluid prop- 
erties across the channel. The magnetic Reynolds number is assumed to be small com- 
pared with unity and thus the induced magnetic field can be neglected. 

SYMBOLS 

The mksA system of units is used herein, 
a quantity defined by equation (17) 

b quantity defined by equation (18) 

B magnetic induction 

C function defined by equation (42) 

d electrode width 
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insulator width 


electric field 

function defined by equation (23) 
function defined by equation (22) 
function defined by equation (24) 

unit vectors in the x-, y-, and z-direction, respectively 

electric current density 

accelerator channel height 

summation indices 

proton charge 

electrode plus insulator width 
magnitude of velocity 
velocity 

distance along X-, Y-, and Z-axes, respectively 
Hall parameter 
current flux function 
efficiency 

convergence parameter 
quantity defined by equation (29) 


scalar conductivity 



r 


mean collision time 


<p electric potential 

<p a constant defined by equation (29) 

(Pfo function defined by equation (29) 

co cyclotron frequency, ^ 

Subscripts: 
e electron 

eff effective 

i ion 

m,n summation indices 

o initial 

T total 

x,y,z x-, y-, and z-components 

Mathematical notations: 

— vector 

unit vector 

| | absolute value 

BASIC EQUATIONS AND ANALYSIS 

For the case of negligible electron pressure gradients and no time dependence, 
the generalized Ohm's law is given by (refs. 1 and 12) 
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( 1 ) 


>ao (E + VXB)-^TxB + ^TxB)xB 


B2 


The last term in equation (1) is the ion slip term and it appears because the ions do not 
transfer their field-acquired momentum with complete efficiency to the neutral particles 
Thus, ion slippage generates heat due to friction and causes, in effect, a decrease in the 
electrical conductivity of the plasma. 

The geometry of the accelerator channel is shown in figure 1. The flow velocity 
of the medium is in the positive x-direction and the applied magnetic field points in the 
negative z-direction; that is, V = Ui and B = Bk where B = B z . It is also assumed 
that all quantities vary only in the x,y plane or -J^ = 0. Thus, equation (1) can be 
written in component form as 

jx , = CT o^x “ 0 e jy - /3 e 0ii x (2) 

jy = °o (®y ” UB) + Peix ~ Pe^Oy ^ 


where /3 e = o> e r e and and o 0 is assumed to be constant not only in the 

axial direction but also across the channel. 


Equations (2) and (3) can be solved for j and j to give 






/3| + (1 + PePi) 


“2[E X (1 + Mi)- (% - UB)/3e 


(4) 


V " ■% , ( r° + ( £ y - DB )( 1 + MO] <«) 

Since the effect of flow velocity gradients and magnetic field gradients along the 
channel is to be investigated, the x-dependence of U and B must be included. Also, 
since (3 e and /3^ are proportional to B, these quantities will also vary with x. As 
has been discussed by Yeh and Sutton (ref. 14), a variation of the velocity across the 
channel does not affect the distribution of E and J (except for a constant factor). 

The problem under investigation is time independent so that the Maxwell equation 
for the conservation of charge becomes 


Let 


- 8jy 

v -S-lf*TiF' 0 

(6) 

' g2 + (J + ^2 

(7) 
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then by using equations (4) and (5) for j x and j y , respectively, equation (6) becomes 


[e x (1 + Mi) - (% - «B)0 e ]£(ln cr eff ) + + ^)(1 + Mi) 


( 8 ) 


where /3 e , /3j, and a gff are assumed not to vary in the y-direction. In obtaining equa- 
tion (8), the Maxwell equation V x E = 0 was used to eliminate two terms. A potential 
( p(x,y ) is now introduced such that 

E = -Vtp (9) 

Equation (8) then takes the form 



d<p 8 (^e^i) _ ch£ 9 ^e 
+ dx 3x 3y 8x 


- ii^e) 


+UB ) 


— In a , f = 0 
9x eff 


( 10 ) 


For the case of constant U, B, and /3, equation (10) becomes the Laplace equation 

3x^ 3y^ 


(11) 


All the previous analyses of the potential and current distribution in MHD accelerators 
did not solve the equation for the electric potential but solved an equation for the current 
flux function y(x,y). Thus if E x and Ey are eliminated from equations (2) and (3) by 
the use of V x E = 0, the following equation results: 



a 

+ V " j x Ifc" + 


CT o 4 (UB) = ° 


( 12 ) 


where V • j = 0 
tion y by 


has been used to eliminate two terms. 


^x 


_ 

~ 9y 


Defining the current flux func- 

(13) 


and 



( 14 ) 
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equation (12) can be written as 




3 y 3 


3 r 


9 0 c 




9(UB) _ Q 
3x 


(15) 


where /3 g and were assumed to be functions only of x. As before, if U, B, and 
|3 are constant, the equation to be solved reduces to the Laplace equation 


8 2 y + 8 2 y 

8x 2 8y 2 


= 0 


(16) 


In the present analysis equation (10) was used to obtain the potential cp. Once q> is 
known, 7 and E can be obtained by using equations (2), (3), and (9). 


The equations for the potential distribution in the accelerator are solved for a 
small section of the accelerator at a time. Therefore, any variation of U and B along 
the accelerator channel can be approximated by expanding these quantities in a Taylor 
series and keeping only terms linear in x. Thus, the variation of U and B becomes 


U = U Q (1 + ax) 

and 

B = B 0 (l + bx) 


(17) 

(18) 


Since and co e are proportional to B, these two quantities vary as 1 + bx. Equa- 
tions (4) and (5) can then be written as 


■'x a o 


E x 1 + /3 e /3^(l + bx) 2 ] - E y /3 e (l + bx) + U 0 B 0 |3 e (l + ax)(l + bx)' 
[l + /3 e /3 i (l + bx) 2 ] + [/3 e (l + bx)] 


(19) 


and 


E x 0 e (l + bx) + E y 1 + ]3 e i3 i (l + bx) 2 ] - U 0 B 0 [l + /3 e /3 i (l + bx) 2 ](l + ax)(l + bx) 

i — : : 


J y 


2 2 
[l + /3 e /S i (l + bx) 2 ] + [/3 e (l + bx)] 


( 20 ) 


where and ^ are the values corresponding to B 0 . Equation (10) then takes the 
form 



( 21 ) 
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where 


g(x) = 


2b/3 e j3^(l + bx) 2b(l + bx)|^2^ e ^ jl + £ e /3j(l + bx) J + /3 e 
1 + ^ e %(l + bx) 2 |i + ^(1 + bx) 2 ] 2 + |8 2 ( 1 + bx) 2 


f(x) = 


b/3g p(l + bx)j2i3ej3i[l + /3 e ^(l + bx) 2 ] + /3J 

1 + ft^U + »x) 2 ; + M . (1 + bx) 2] 2 + ^ + bx) 2 


/ \ U o B o^e(l + bx ) / 
h(x) = ^ a + 3abx + 2b 

1 + j3 e /%(l + bx) 2 l 


2b(l + ax)(l + bx) 2 jj2/3 e /3j|] + @^(1 + bx) 2 ] + /3 2 
[l + i3 e ^(l + bx) 2 ] + /3 2 (1 + bx) 2 


To obtain the potential distribution in the channel, equation (21) was solved numeri- 
cally with the appropriate boundary conditions. At the surface of the conducting elec- 
trodes, the potential gradient is taken to be zero; thus, 


E x = 0 


<p = Constant 

over the electrode surface. From equations (4) and (5), one then finds that 

3x _ -/M 1 + bx ) 

3y 1 + 0 e 0j(l + bx) 2 

over the electrode surface. The current into an insulator is zero; thus 

jy = ° 


' — 2 I* + foT + ^BoC 1 + ax)(l + bx) = 0 

i + jSg/^ci + bx) 2 
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over the surface of an insulator. One also finds that 


E„ -J3 P (1 + bx) 

x _ ' /ggx 

Ey - UB 1 + /3 e /%(l + bx) 2 

over the insulator surface. 

In the present analysis the calculations are performed for the case where each pair 
of electrodes has a separate power supply. The cathode and corresponding anode can be 
shifted any amount along the channel but the current out of dhe electrode equals the cur- 
rent into the other electrode of the pair. This additional boundary condition is then used 
to fix the potential along the two remaining boundaries as indicated in figure 2 by means 
of an iterative method. For example, if a = b = 0, then the velocity and magnetic field 
gradients along the channel are zero and cp(x, y) - cp(x - 2s, y) is a constant, where s 
equals the sum of the length of one insulator plus one electrode. A section of the genera- 
tor of length 2s was chosen for all the numerical solutions of equation (21). If a or 
b is not equal to zero then cp(x - 2s, y) - <p(x,y) = ^(y)- The expression for ^(y) that 
was used in most of the calculations is of the form 

<P h (y) = <p a + py (29) 

where <p a and p are constants that must be chosen to satisfy the boundary conditions. 

It is of interest to determine the fraction of the power input that is used to accel- 
erate the fluid. Let e be the ratio of the work done per unit time by the Lorentz force 
in the flow direction to the power input; thus, 


JJ V • (I x b)<±x dy 
JJ J . E dx dy 


(30) 


where e as defined by equation (30) is taken to be the efficiency of the accelerator. For 
the present problem, this ratio becomes 


e 



dx dy 


jJ(3x E x + jy E y) dx d y 


(31) 


In the numerical calculations, equation (21) was used and only a z-component of the 
magnetic field was considered. Since it is assumed that B is obtained by external 
excitation, the relation 

V x B = 0 (32) 
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or 


9B Z _ 9B X 
9x 9z 


( 33 ) 


must be satisfied. For the variation of B z given in equation (18), one then obtains for 
B x the variation 

B x = bB 0 z (34) 

where the constant of integration is taken to be zero. The total magnetic field is then 
given by 

B = ibB 0 z + kB 0 (l + bx) (35) 

Instead of equations (2) and (3) the following equations for the components of equation (1) 
are obtained: 

/3„ ( 2 \ 

lx = a o^x - "g" Bziy + g2 \ x®z3z “ ®z^xj 

p e . 

iy = °o(®y - UB Z ) - "g'^B x j z - B z j x J- (37) 


h = °o E z + § B x j y + ^(B x B z j x - 


(38) 


The solution of this system of equations is too difficult by a numerical method. To 
simplify the equations, it can be assumed that j = 0 since no current can flow through 
the insulated side walls. From equation (38), such a condition requires a z-component 
of the electric field given by 


E z = 


ffeBx . ffeftBxBz . 
B°o 3y " B 2 a Q ]x 


(39) 


However, near the electrodes E z is shorted so that, in general, the z-component of the 
current density cannot be zero. Thus, in addition to the concentration of current at the 
upstream edge of the cathode and at the downstream edge of the anode obtained later in 
the section on results, the x- component of the magnetic field will now cause additional 
current concentration. Depending on whether b is positive or negative, the current 
will either concentrate at the center of the electrode edge or at the two ends of the elec- 
trode edge. If b or z is small, then the effect of B x can be assumed to be small 
and equations (21) to (24) describe the system. In this manner some information about 
the effect of the axial magnetic field gradient can be obtained even though some simpli- 
fying assumptions had to be made. For the case of an axial velocity gradient, no such 
difficulties exist and the system will remain so that it can be treated as two dimensional. 
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NUMERICAL SOLUTION 


Equation (21) can be written in the following finite difference form: 


^m+1. 


(if . g(x)l . w [(if §m1 _ l0 

,n[_\Axy + 2 Axj + ^m-l,n[\Ax/ 2 AxJ ^m,n 


2 , 2 
(Ax ) 2 (Ay)' 


+ V,n+1 


Ayy 2Ay 


+ V 


m,n - 1 


JL\ _ l&L = h(x) 

Ayy 2Ay W 


( 40 ) 


where Ax and Ay represent the separation of the mesh points used in the numerical 
solution. The subscripts m,n refer to the mesh points; n = 1 represents a line of 
points along the bottom row of electrodes and insulators; n = N represents the points 
along the top row of electrodes and insulators of the accelerator channel; m = 1 is the 
leftmost column of mesh points; and m increases with increasing x. 

At points along the insulator j y = 0 and equation (27) can be written in finite dif- 
ference form as 


2Ax (Vn + l ' V,n-l) + * Vl,») " ‘ 4Ax A y U o B o» + + bx > < 41 > 

where 

|3„(1 + bx) , „ 

C(x) = — -37 (42) 


1 + /3 e /3i(l + bx)' 


Equation (41) is used to eliminate the <p or the cp m n+ ^ term which appears in 
equation (40) whenever any points fall on the insulators. Thus for any points along the 
bottom insulators, n = 1 and 


* s - - - iiH-ufe) 2 - ig+m* - $ 


2 2 
(Ax ) 2 (Ay ) 2 


= h(x) + U 0 B 0 (1 + ax)(l + bx) 


(43) 


must be satisfied. Also, along the top insulators n = N and the relation 


11 



‘ p m+l,N{^I + |i|- [ f < x> + ^} + 'Ita + Sf [ 1<x) + 


- v»^! + <^js] * Vn-i[^] - h « * U 0 B 0 (1 + «*)<! + bx) 


f(x) +-£- 
Ay 


( 44 ) 


must be satisfied. The boundary condition along the electrodes is cp m ^ = Constant and 
cp m N = Constant. Along the two remaining boundaries, the condition given by equa- 
tion (29) 


^n " ^M,n 


cp^-PY 


(45) 


is used. The quantities cp„ and p in equation (45) must be chosen so that the current 

a 

from an anode goes into the corresponding cathode. The correct values of p and cp 

ci 

are determined by integrating the current density along a line between opposing insula- 
tors and satisfying the condition 


S 


j x dy = pI T 


(46) 


where I<p is the total current per unit width of the channel through one electrode and 
p = 0, ±1, ±2 depending on the amount of electrode stagger. If U and B are constant, 
then p = 0 and only cp a needs to be determined. 

The system of simultaneous equations given by equations (40), (43), and (44) were 
solved on the IBM 7094 data processing system by an iteration method. The form of 
equation (40) actually put on the computer is 


q Q r+1 = cp r + v(r)l 


iffJ-f + (-L)* 

VAx/ \Ay/_ 


, w - [?JLf + hM] . ^ r +l 

y m+l,n[vAx/ + 2AxJ 


L 

l,n _\Ax/ 




2Ax 


+ Vi 






Ay/ 2Ay 


+ cp r+1 1 


iAyy 2Ay 


- h(x) > - (p 


m,n 


(47) 


Equations (43) and (44) can be done similarly. The superscript r refers to the rth 

iteration and v is a convergence factor which was used to speed convergence and could 

be varied during the solution. For points along the insulators, an equation of the form 

of equation (47) was initially tried. However, as found by previous workers (refs. 8 

and 15), this form of the difference equation caused the numerical solution to blow up. 

r+1 


For the points along the insulator, any change in cp 


from the previous iteration was 
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then divided by 1 + j8 e and fairly rapid convergence was obtained. Another difficulty 
which was due to the large potential gradients near one edge of the electrodes was 
encountered. Again this difficulty was primarily for points on the insulator. For the 
first few mesh points near the electrode edge, there was an oscillation (in space) of the 
potential such that every other cp m j or <p mjN deviated from the "correct" <p. This 
problem could be reduced but not eliminated by going to a finer mesh size. The difficulty 
was overcome by checking for the oscillations near the electrode edge during the solu- 
tion. If the oscillations in cp became appreciable, an averaging of the potential was 
performed to eliminate the oscillations and the solution was continued until the system 
converged. 

The total power input to the accelerator section was computed by summing 


where j x and jy 


and 


M N 


I 1 [W'-k. + ( i y)m,n( E y)m,n] A * 
m=l n=l 

are given by equations (19) and (20) and 

(48) 

^^x)m,n 2Ax(^ni-l,n ^m+l,n ; 

) 

(49) 

(^Vjmjn = 2Ay (^m,n-l " ^m,n+l) 


(50) 


RESULTS 


The parameters used for most of the calculations were those that are pertinent 
for the 1-inch- (2.54-cm) square accelerator presently in operation at the Langley 
Research Center. These parameters are given in reference 11. To determine the 
effect of electrode stagger, the following parameters were used: a Q = 500 mho/ m, 

B 0 = -1.14 webers/m^, p e = -16, and U = 400 m/sec in the midsection of the accel- 
erator channel. The channel is described by the following geometrical parameters: 
l = 2.54 cm, d = 0.72 cm, and e = 0.41 cm. (See fig. 2.) The applied potential between 
corresponding electrodes was 200 volts. In order to show the effects of ion slip, the 
calculations were made for both j3j = 0 and /3j = -0.04 with a = b = 0. The results 
are shown in figures 3 to 7. The shift position has been normalized to s, the length of 
one insulator width plus one electrode width. An increase in the shift position corre- 
sponds to moving the anode in the flow direction. In figure 3 the calculated Hall field 
is plotted as a function of electrode shift position. The results show that the Hall field 
peaks between the 0 and - 1 shift position and that ion slip appreciably reduces the Hall 
field. Figure 4 shows the effect of electrode stagger on the current per electrode. As 
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expected, ion slip reduces the current. In contrast to the Hall field, the current peaks 
at the 0 shift position. Similarly, figure 5 shows that the calculated internal resistance 

^applied + BUZ 

I T l 


between a cathode and the corresponding anode is a minimum at the 


0 shift position. The equation for the internal resistance given by Witalis (eq. (42) in 
ref. 5) was also used to calculate the resistance for zero shift. The results were found 
to be identical to those given in figure 5 for zero shift. The variation as a function of 
electrode stagger of the power per unit volume used in accelerating the fluid in the axial 
direction is shown in figure 6. The acceleration peaks at the zero shift position and at 
that point ion slip reduces the work done by about 7 percent. Figure 7 shows the cal- 
culated effect of electrode stagger on efficiency. It can be seen that the efficiency peaks 
for zero shift. Ion slip effects reduce the efficiency by only about 2 percent. Thus, the 
calculations show that shifting the anodes relative to the corresponding cathodes will 
reduce the efficiency and acceleration. Thus, on the basis of the present calculations, the 
shifting of the electrodes which has been proposed to increase acceleration (ref. 16) can- 
not be recommended. Carter et al. (ref. 16) have experimentally investigated the effect 
of shifting the electrodes. Their results are shown in figures 8 to 10. The experi- 
mentally obtained Hall field and current per electrode decreases much more rapidly with 
shift position than the calculated values. This result is probably due to the effects of 
the boundary layer near the electrodes and insulators which may cause much stronger 
variation of the total internal resistance with electrode shift than that calculated. Fig- 
ure 8 shows a peak in the Hall field between the 0 and - 1 shift position which is in agree- 
ment with the calculated results. The scatter of the data in figures 9 and 10 is very large 
and, as a result, it is difficult to compare these experimental results with the calculated 
values given in figures 4 and 6. 


The equipotential and current flux lines in the accelerator section are also of 
interest. In figure 11 the equipotential lines are shown for the different shift positions 
for the case with /3j = -0.04. Figure 12 shows the lines for constant current flux, again 
for /3j = -0.04. The difference between the adjacent equipotential lines A cp and adja- 
cent current flux lines Ay are as given in the figures. The curves were plotted by 
interpolating between the calculated potential at the different mesh points. Also, since 
the lines of constant current flux were obtained from the potential by using the approxi- 
mate formulas for E in equations (4) and (5), one cannot expect very high accuracy of 
the slope of these curves in the region of high current density concentration near the 
edge of the electrodes. All the equipotential and constant current flux curves were 
obtained by starting in the lower left-hand corner of the accelerator section where 
y = y 0 and cp = cp Q , and by then finding curves of constant cp and y given by 
cp = cp Q ± n A cp and y = y Q ± n Ay where n = 1, 2, 3, . . . . 


I 
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An attempt was made to calculate a case which matched the experimentally 
obtained data given in reference 17. Since the only accurate data available are the aver- 
age current per electrode and the Hall field, it is not possible to obtain a unique match. 

For example, for given arbitrary values of /3 e and /3j, the applied field Ey can be 
chosen to give the desired Hall field and a Q can then be chosen to give the desired 
current per electrode. It was therefore assumed that /3 e = -16 and = -0.04 are 
the correct values for the Hall parameter and the ion slip term, and the applied potential 
and cr 0 were chosen to give the experimentally obtained Hall field and current per elec- 
trode as given in reference 17. The magnetic field was -1.14 weber s/m^ and the flow 
velocity was assumed to be 2000 m/ sec at the channel entrance and 6000 m/sec at the 
exit and to vary linearly in between. Then o Q was required to be 355 mho/m. In fig- 
ures 13 to 15 the results of the calculations are compared with the experimental data 
given in reference 17. Figure 13 shows the measured current per electrode and the 
current used in the calculations. In figure 14 which shows the calculated and the mea- 
sured potential between electrodes, the difference between the experimental points and 
the calculated curve should represent the sheath drop near the electrodes and insulators. 
The figure indicates that either o Q increases or the sheath drop decreases near the 
end of the channel. In figure 15 the calculated electrode potentials and their variation 
along the accelerator (Hall potential gradient) are compared with the experimental data 
given in reference 17. The agreement is very good. 

The effect of /3 e on the electrical parameters of the accelerator has also been cal- 
culated. For these calculations was taken to be zero and the voltage applied between 
an anode and corresponding cathode was 200 V. First the results for the case with 
B = Constant = -1.14 webers/m^ and cr = Constant ■ 500 mho/m are given. Since B 
and a are constant, /3 e will vary as j3 * Constant/n where n is the electron density. 
Figures 16, 17, and 18 show the variation with j3 e of the Hall field, the current per 
electrode, and the efficiency. The variation of the internal resistance is shown in fig- 
ure 19. Figure 16 shows that for -j3 e = 16 the Hall field is somewhat lower than that 
shown in figure 15. The reason for this difference is, of course, that the applied poten- 
tial (or Ey) is lower for the case shown in figure 15. As shown in figure 17 for constant 
applied voltage, the current per electrode decreases rapidly with increasing Hall param- 
eter. Also shown in figure 17 is the internal resistance as calculated from equation (42) 
of reference 5. It can be seen that the two methods give identical results except for 
-j3 e larger than about 20. 

When /3 e varies such that j3 e = 14B and a Q = Constant = 500, the results are 
as shown in figures 20 to 23. The efficiency increases almost linearly with -/3 e but the 
work done in accelerating the fluid decreases rapidly. For sufficiently large -/3 e , the 
accelerator would turn into a generator but this transformation occurs only because a 
constant applied voltage was assumed, and as B increases, so does the 
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U x B . I voltage. For large enough B, this voltage exceeds the applied potential. 

The internal resistance for this case is the same as that shown in figure 19. 

In figure 24 the variation of the internal resistance as a function of the ratio of 
electrode width to s is shown. The solid line is obtained from the equation for the 
internal resistance given by Witalis (eq. (42) of ref. 5) and the circles are the values 
computed by the present numerical method. For the results shown in figure 24, 

B = -1.14, a Q = 500, fa = -0.04, s = 0.113 cm, and l = 2.54 cm. Again, there is 
excellent agreement for the two methods. 

In an accelerator the velocity along the channel increases and the quantity a in 
equation (17) is not zero. In figures 25 and 26 the internal resistance and the efficiency 
are shown as functions of a. The parameters used are the same as for the previous 
case with /3 e = -16 and d = 0.72 cm. The results shown in figures 25 and 26 show that 
the internal resistance increases and the efficiency increases with increasing velocity 
gradients, but both change only slightly over a rather large range of gradient. For 
a = 3.33 (the operating conditions of the 1-inch- (2.54-cm) square accelerator at the 
Langley Research Center), the equipotential and constant current flux plots are nearly 
identical to those for a = 0. However, for larger a the equipotential lines and 
constant current flux lines begin to curve, as shown in figure 27. The Hall field along 
the center line of the accelerator remained constant for increasing a. However, p in 
equation (29) had to be increased to satisfy the boundary conditions. Figure 28 shows 
py N as a function of a, and indicates how the voltage applied between electrodes must 
be increased along the accelerator channel. 

The effect of magnetic field gradients with the approximations stated earlier has 
also been investigated. The results are shown in figures 29 to 33. Figure 29 shows that 
the internal resistance increases slightly with increasing b. As can be seen from fig- 
ure 30, the efficiency increases as b increases. The change in the equipotential and 
constant current flux plots is shown in figures 31 and 32. As the magnetic field gradients 
become large, either negatively or positively, the current density and electric field 
become increasingly more nonuniform. The effect of b on the applied voltage is indi- 
cated in figure 33 which shows py^ as a function of b. 

For applied magnetic field gradients of the order of the self-induced field, the 
numerical solution is almost identical to that for a zero magnetic field gradient. This 
result justifies the neglect of the self- induced field. For increasing velocity and mag- 
netic field gradients, it becomes more difficult to obtain convergence of the numerical 
solution. The numerical solution could probably be improved by using equation (28) 
instead of equation (27) as the boundary condition over the insulator surface. 

All results presented so far were for the 2. 54 -cm- square accelerator. Calcula- 
tions were also performed to determine the effect of electrode stagger for a larger 
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accelerator with a different ratio of if s. For this accelerator, l = 6.35 cm, 
e = 0.476 cm, and d = 0.795 cm. The other parameters used for the calculation were: 

B = 1.25 webers/m 2 , /3 e = -10, a Q = 500 mho/m, U = 7000 m/sec, and a potential dif- 
ference of 640 volts between corresponding electrodes. The results were found to be 
very similar to those obtained for the smaller accelerator; consequently, only the effect 
on the Hall field and the internal resistance are shown (figs. 34 and 35). The internal 
resistance varies very slowly with shift position. 

CONCLUDING REMARKS 

The calculated effects of electrode stagger on the electrical parameters of an 
accelerator were found to be in qualitative agreement with experimental results. It was 
found that shifting the electrodes will reduce efficiency and acceleration. The more 
rapid decrease of the parameters found experimentally may be due to the neglect of the 
variation of the fluid parameters in the vertical (y) direction. The numerically obtained 
variations of the internal resistance of the accelerator with Hall parameter and ratio of 
electrode to electrode plus insulator width (d/s) for zero velocity and magnetic field 
gradients were in agreement with the results obtained by Witalis by means of a conformal 
method. 

With increasing axial velocity and magnetic field gradients, the current density and 
electric field distribution become increasingly nonuniform in the channel center. How- 
ever, the variation of the electrical parameters with velocity or magnetic field gradients 
obtained by the present method indicates that there are no serious detrimental effects on 
accelerator operation for moderate gradients. Whether problems arise because of addi- 
tional current concentration due to the appearance of the magnetic induction (B x ) com- 
ponent for axial magnetic field gradients must still be analyzed. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 9, 1967, 

129-02-01-01-23. 
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Anode shift position relative to cathode, s Anode shift position relative to cathode, s 

Figure 5.- Effect of electrode stagger on the internal resistance between electrodes. Figure 6.- Effect of electrode stagger on work done to accelerate the flow. 
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Figure 7.- Effect of electrode stagger on accelerator efficiency. 
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Figure 8.- Experimentally obtained variation of Hall field with electrode stagger taken from reference 16. 
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Figure 13.- Comparison of experimentally obtained electrode current with current used in calculations. 
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Figure 14.- Comparison of calculated potential applied between electrodes with experimental values. 
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Figure 15.- Comparison of calculated electrode potential with experimental data from reference 17. Figure 16.- Hall field as a function of the Hall parameter with B and o constant. 
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Figure 25.- Internal resistance as a function of velocity gradient. 
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Figure 26.- Variation of efficiency with velocity gradient. 
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Equipotential lines j Equipotential lines; 

a =6; Acp = 10 V a =10; Acp=10V 


Lines of constant current flux* Lines of constant current flux 

a = 6; a y = 200 A/m a = 10; a y - 200 A/m 




Figure 27.- Variation of equipotential lines and current flux lines with a. The top electrodes are the anodes and the flow 

is from left to right. 







Figure 28.- Variation of pY^ with a. 
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Figure 29.- Internal resistance as a function of b. 
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Figure 30.- Efficiency as a function of magnetic field gradient. 
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Figure 32.- Variation of constant 
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Figure 33.- Variation of pY N with b. 
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Figure 34.- Variation of Hail field with electrode stagger for the 6.35-cm-square accelerator. 
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Figure 35.- Variation of internal resistance with electrode stagger for the 6.35-cm-square accelerator. 
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